Paradoxical transitions to instabilities in hydromagnetic Couette- Taylor flows 
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By methods of modern spectral analysis, we rigorously find distributions of eigenvalues of lin- 
earized operators associated with an ideal hydromagnetic Couette- Taylor flow. Transition to insta- 
bility in the limit of vanishing magnetic field has a discontinuous change compared to the Rayleigh 
stability criterion for hydrodynamical flows, which is known as the Velikhov-Chandrasekhar paradox. 
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Instabilities of Couette- Taylor (CT) flow between two 
rotating cylinders are at the cornerstone of the last- 
century hydrodynamics In 1917 Rayleigh found a 
necessary and sufficient condition for centrifugal insta- 
bility of CT-flow of an ideal fluid between cylinders of 
infinite length with respect to axisymmetric perturba- 
tions [2]. Taylor extended Rayleigh's result to viscous 
CT-flow and computed seminal linear stability diagrams 
that perfectly agreed with the experiment at moderate 
angular velocities [3|. 

Despite the theoretical and experimental studies of the 
Couette- Taylor flow are more than a century long, re- 
cent decade had seen a true renaissance of this classi- 
cal subject caused by the increased demands of the ac- 
tively developing laboratory experiments with liquid met- 
als that rotate in an external magnetic field [4]. The 
prevalence of resistive dissipation over viscous dissipa- 
tion in liquid metals dictates unprecedentedly high val- 
ues of the Reynolds number (Re ~ 10 6 ) at the threshold 
of the magnetorotational instability (MRI) of hydrody- 
namically stable quasi-Keplerian flows that currently is 
considered as the most probable trigger of turbulence in 
astrophysical accretion discs Difficulties in keeping 
hydrodynamical CT-flows laminar at such high speeds, 
put the laboratory detection of MRI at the edge of mod- 
ern technical capabilities. 

Is the existing theory of MRI well-prepared to face 
these promising experimental opportunities? No matter 
how paradoxical it may sound, the answer is: Not yet. 

Indeed, already the discoverers of MRI, Velikhov Q 
and Chandrasekhar pj], pointed out a counter- intuitive 
phenomenon. In case of an ideal non-resistive flow, which 
we consider in this Letter, boundaries of the region of the 
magnetorotational instability are misplaced compared to 
the Rayleigh boundaries of the region of the centrifugal 
instability and do not converge to those in the limit of 
negligibly small axial magnetic field. In presence of dis- 
sipation the convergence is possible 0. 

The existing attempts of the physical explanation of 
the Velikhov-Chandrasekhar paradox [9j involve Alfven's 
theorem that 'attaches' magnetic field lines to the fluid 
of zero electrical resistivity, independent of the strength 



of the magnetic field, which implies conservation of the 
angular velocity (Velikhov-Chandrasekhar) rather than 
the angular momentum (Rayleigh). However, the weak 
point of this argument is that the actual boundary of 
MRI does depend on the magnetic field strength even in 
the case of ideal MHD and tends to that of solid body 
rotation only when the field is vanishing. This indicates 
that the roots of the paradox are hidden deeper. 

Recently, this intriguing effect was reconsidered in the 
full viscous and resistive setting by a local WKB approx- 
imation It was found that the threshold surface 
of MRI in the space of resistive frequency, Alfven fre- 
quency and Rossby number possesses a structurally sta- 
ble singularity known as the Plucker conoid that persists 
at any level of viscous dissipation. The singular surface 
connects the Rayleigh- and the Velikhov-Chandrasekhar 
thresholds through the continuum of intermediate states 
parameterized by the Lundquist number [Icj |. 

Why does this singularity exist? Our Letter sheds light 
to this question via rigorous inspection of the spectra of 
the boundary eigenvalue problems associated with the 
ideal hydrodynamic and hydromagnetic CT-flows. Rig- 
orous spectral results are illustrated by MATLAB com- 
putations of eigenvalues of the linearized operators. 

If u is the velocity field, b is the magnetic field, and 
cylindrical coordinates (r, 6, z) are used, the basic CT- 
flow between cylinders of radii R\ and R2 , Ri < R2 , is 
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where &o is arbitrary and (a, c) are related uniquely to 
O1.2 = £l(Ri,2) through the viscous limit, 
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In the case of co-rotating cylinders, Q.%,2 > 0, the 
Rayleigh boundary corresponds to a = 0, whereas the 
Velikhov-Chandrasekhar boundary is c = 0. 

The summary of our results is as follows. 

(I) In the case of no magnetic field (bo — 0), co-rotating 
cylinders (fii,2 > 0), and an ideal fluid, we prove that 
the linearized stability problem has a countable set of 
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neutrally stable pairs of (purely imaginary) eigenvalues 
for a > and a set of unstable pairs of (purely real) 
eigenvalues for a < 0, all accumulating to zero. At a = 0, 
all pairs of eigenvalues merge together at zero. 

(II) Under the same conditions but for counter-rotating 
cylinders with Sli < and 0,2 > 0, we show that there 
exist two sets of eigenvalue pairs: one set contains real 
eigenvalues and the other set contains purely imaginary 
eigenvalues. The unstable real eigenvalues converge to 
the zero accumulation point when fix — > for fixed 
Q.2 > (where a > 0), whereas the stable imaginary 
eigenvalues persist across Oi = 0. 

(III) For any magnetic field (&o 7^ 0), co-rotating cylin- 
ders (fii ( 2 > 0), and an ideal non-resistive hydromagnetic 
flow, we prove that there exist two sets of eigenvalue pairs 
and both sets contain only purely imaginary eigenvalues 
for < il± < f^. One set remains purely imaginary 
for fii > il 2 but the other set transforms to the set of 
real eigenvalues along a countable sequence of curves, 
which are located for il\ > H2 and approach the diago- 
nal line ili — ^2 (c = 0) in the limit 60 0. One pair 
of purely imaginary eigenvalues below the corresponding 
curve transforms into a pair of unstable real eigenvalues 
above the curve. No eigenvalues pass through the ori- 
gin of the complex plane in the neighborhood of the line 
a — 0, even if bo is close to zero. 

(IV) Under the same conditions but for counter- 
rotating cylinders with f2i < and O2 > 0, we show 
the existence of four sets of eigenvalue pairs, which are 
either purely imaginary or real. The unstable eigenval- 
ues bifurcate again along a countable sequence of curves, 
which are located for f2i < and approach Q± = in 
the limit bo — > 0. The purely imaginary pair of eigen- 
values above the curve turns into a purely real pair of 
eigenvalues below the curve. 

Although the results (I) and (II) partially reproduce 
the conclusions of Synge , the existence of zero eigen- 
values of infinite multiplicity at the Rayleigh threshold is 
emphasized here for the first time. Similar coalescence of 
all eigenvalues at the zero value happens also in the Bosc- 
Hubbard dimer [12j. Results (III) and (IV) are new to 
the best of our knowledge. Numerical evidences of these 
results can be found in [13| ]. 

The rest of our paper is devoted to the proofs of the 
above results and their numerical illustrations. We take 
the equations for an ideal hydromagnetic fluid 0, 



u t + (u • V)u = -V (p ■ 
b( = Vx(uxb), 
V-u = 0, V-b = 0, 
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(3) 



where p is the pressure term determined from the incom- 
pressibility condition V • u = 0. We linearize (J3j> at the 
basic flow (flj and use the standard separation of vari- 
ables for symmetric (^-independent) perturbations, 



where 7 is the growth rate of perturbations in time and 
fc 6 1 is the Fourier wave number with respect to the 
cylindrical coordinate z. Performing routine calculations 
8] , we find the system of four coupled equations for com- 
ponents of U and B in the directions of e r and eg (de- 
noted by U r , Ug, B r , and Bg), 



ikb (k 2 + L)B,, 
ikboBg — 2aU r 
ikboU r = jB r , 
ikb Ug - ^B r = 1 B, 
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7 (fc 2 + L)U r , 



(5) 



u = u + U(r)e 
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b = b + B(r)e 
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where L = —d 2 — ^d r + \ is the Bessel operator, which 
is strictly positive and self-adjoint with respect to the 
weighted inner product (/, g) — J^ 2 rf(r)g(r)dr. We 
note that z-componcnts of U and B, as well as the pres- 
sure term, have been eliminated from the system of equa- 
tions ([5]) under the condition k ^ 0. 

For hydrodynamic instabilities of the CT-flow, we set 
bo = 0, which yields uniquely B r = Bg = 0, 2aU r +jUg = 
0, and a closed linear eigenvalue problem, 

7 2 (/c 2 + L)U r = -Ak 2 aVL{r)U r , R t < r < R 2 , (6) 

subject to the Dirichlet boundary conditions at the inner 
and outer cylinders U r (R\) = U r (R2) = 0. 

The operator L is an unbounded strictly positive op- 
erator with a purely discrete spectrum of positive eigen- 
values {/i„}„gt( that diverge to infinity according to the 
distribution p, n ~ n 2 as n — > 00. Inverting this operator 
for any k £ IR and defining a new eigenfunction ^ by 
U r = (k 2 + L)~ 1/2 <I>, we rewrite © in the form, 

7 2 * = -aT*, T = 4k 2 (k 2 +L)- 1/2 n(k 2 +Ly 1/2 , (7) 

where the self-adjoint compact operator T has eigen- 
values {— 7 2 /a}„ S N that accumulate to zero with 7„ = 
C(n _1 ) as n — > 00. 

If fii,Q 2 > 0, then fi(r) > for all r G [Ri,R 2 ] and 
T is a compact positive operator. Hence, all 7 2 < if 
a > and all 7 2 > if a < 0. The condition a = 
(^2^2 = ^l-R?) is the Rayleigh boundary, at which all 
eigenvalues are at 7 = 0. The proof of (I) is complete. 

If Oi < and 2 > 0, then a > but f2 is 
sign-indefinite on [i?i,i?2]- Since T is a compact sign- 
indefinite operator, it has two sequences of eigenvalues 
accumulating to zero: one sequence has 7 2 < and the 
other one has 7 2 > 0. This completes the proof of (II). 

Figure [lja) gives numerical approximations of the five 
largest and five smallest squared eigenvalues 7 2 as func- 
tions of the parameter VL\ for fixed values of = 1, 
R\ = 1, i?2 = 2, and k = 1. The dotted line shows the 
accumulation point 7 = for the sequences of eigenval- 
ues. For ill > 0, the five smallest eigenvalues are not 
distinguished from the zero accumulation point. 

For hydromagnetic instabilities, we express B r , Bg, 
and Ug from the system of linearized equations ([S]) and 
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(a) 





(b) 



FIG. 1: Squared eigenvalues j 2 of the problem (0 versus 
fii for 0,2 = 1, Ri = 1, R 2 = 2, and k = 1. (a) & = 0: 
all squared eigenvalues 7 coalesce to zero at the Rayleigh 
line fii-Ri = Q.2R2, whereas positive squared eigenvalues for 
fii < merge to zero at fii = 0. (b) bo = 0.4: the squared 
eigenvalues of the problem ([HJ> change stability above the 
Velikhov-Chandrasekhar line fii — Q2 and below fii = 0. 



find a closed linear eigenvalue problem, 



aY U, 



{ 1 2 + k 2 b 2 ) 2 (k 2 + L)U r =Ak 2 n(r)[^ 

(8) 

subject to the same Dirichlet boundary conditions at 
r = Ri t 2- If bo = and 7 ^ 0, system (|5J) reduces to 
©, however, it is a bi-quadratic eigenvalue problem and 
hence has a double set of eigenvalues compared to (JB]). 

Denoting A = 7 2 + fc 2 6g, we rewrite © as the quadratic 
eigenvalue problem, 

A 2 (fc 2 + L)U r + Aak 2 \VL{r)U r = 4fc 4 6 2 ^ 2 (r)J7 r . (9) 

It follows again from the compactness of the operators 
(fc 2 + L) _1 ri and (fc 2 + L^Vi 2 that the spectrum of 
the quadratic eigenvalue problem © is purely discrete. 
Chandrasekhar [7] showed that all eigenvalues A are real. 
We shall prove that these eigenvalues accumulate to zero 
as two countable sets with A„ = C(n _1 ) as n — > 00, one 
set is for positive A and the other set is for negative A. 
The result definitely holds for a = because A 2 becomes 
an eigenvalue of the self-adjoint problem, 

A 2 * = fc 2 ^*, S = Ak 2 {k 2 + L)- 1/2 n 2 (k 2 + L)~ 1/2 , 

(10) 

where S is a compact positive operator. 



To show the same conclusion for a ^ 0, we use a re- 
cently developed technique from [14J and rewrite (j9)) as 
a parameter continuation problem for v = A -1 , 



avVt{r)U r 



(k' z + L)U r + k 2 ble 2 n 2 (r)U r . (11) 



4fc 2 



Here eigenvalues v of (jTTJ) for are continued with 
respect to the real values of e to recover eigenvalues A = 
v^ 1 of ((9]) at the intersections with the diagonal v = e. 

At e = 0, we recover back the hydrodynamical problem 
©. If n 1( Cl 2 > 0, then n(r) > for all r G [Ri,R 2 ] and 
eigenvalues {^Vi(e)}neN at e = are strictly negative if 
a > or strictly positive if a < 0. Moreover, v n (0) ~ n 2 
as n — > 00. Without loss of generality, let us consider 
the case a > 0. Each negative eigenvalue v n {e) is strictly 
increasing for large values of |e| at any point eo, because 



ae - 



de 



€=e 



>0, 



(12) 



where tp n is the eigenfunction for the eigenvalue ^n(e) 
in (fTT|) at e = eo. The right-hand-side of (|12|) is always 
bounded, hence the eigenvalues {fn(e)}neN are contin- 
ued to positive infinity as |e| — > 00. As a result, there 
exist two countable sets of intersections of eigenvalues 
{i/„(e)} ra gN with v = e, one set is for positive A = v~ x 
and the other set is for negative A. Both sets accumulate 
at zero as n — > 00. This completes the proof of (III). 

If ill < and r?2 > 0, then a > but f2 is sign- 
indefinite on [Ri, R 2 ]. In this case, again using the com- 
pact operator T in (JTJ), there exist two sets of eigenval- 
ues {z^(e)}neN of (HO at e = 0: one set {^(0)}„ £ n 
is strictly negative with {flip^ ,fn) < an d the other 
set {i/+(0)} n£ N is strictly positive with (Q<^+,<p+) > 0. 
Because the signs of (Q<p^,(p^) are preserved for small 
e 7^ 0, it follows from the derivative (fT2]) that the eigen- 
values {v~ (e)}„ S N are convex upward for larger values of 
|e| and the eigenvalues -{V+(e)}raeN are concave downward 
for larger values of e. The curves of {z^(e)}„ S N may in- 
tersect but the intersection is safe (i.e., eigenvalues split 
without onset of complex eigenvalues) because the eigen- 
value problem (jlip is self-adjoint for any real e and hence 
multiple eigenvalues are always semi-simple. If the signs 
of {SltPni fn) are preserved along the entire curves, then 
we conclude on the existence of four sets of intersections 
of these eigenvalues with the main diagonal v = e: two 
sets give positive eigenvalues A and the two other sets 
give negative eigenvalues. The conclusion is not affected 
by the fact that (Tl(p^,<p^) may vanish along the curve. 
If this is happened, then (Qip^ , tp^) has at least a simple 
zero due to analyticity in e and hence the derivative (|12[) 
implies that the corresponding curve ^(e) goes to plus 
or minus infinity for finite values of e. This argument 
completes the proof of (IV). 

Figure [ljb) shows numerical approximations of the five 
smallest and five largest squared eigenvalues 7 2 as func- 
tions of J7i for fixed values of il 2 = 1, Ri = 1, R 2 = 2, 
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FIG. 2: Curves of zero eigenvalues (the stability domain is 
located between red curves) for R\ = 1, R 2 = 2, bo = 0.4, 
and k = 1. The curves approach the line fli = SI2 as bo — > 0. 

60 = 0.4, and k = 1. Cascades of instabilities arise for 
Sli > O2 and fii < by subsequent merging of pairs of 
purely imaginary eigenvalues 7 at the origin and splitting 
into pairs of real (unstable) eigenvalues 7. For Qj > 0, 
the two sets of squared eigenvalues accumulate to the 
value 7 2 = — k 2 bg (A = 0), which is shown by the dotted 
line. For Sli < 0, a more complicated behavior is ob- 
served within each set: the squared eigenvalues coalesce 
and split safely, indicating that each set is actually rep- 
resented by two disjoint sets of the squared eigenvalues. 

To study the instability boundaries in ([8]). we substi- 
tute 7 = and regroup terms for ba ^ to obtain 

bl(k 2 + L) Ur = 4(0! - n 2 ) 1 &^u r , (13) 

If O12 > 0, it follows from equation (113)) that there 
exists a countable set of bifurcation curves for Or > f^, 
because L is a positive operator and Cl(r) is strictly pos- 
itive. On the other hand, in the quadrant fli < and 
> 0, there exists another set of bifurcation curves, 
because il is sign-indefinite and L is unbounded. 

To study further the instability boundaries, we notice 
that fi(r) depends on both J7i and f^- Therefore, we 
shall rewrite ([TBf as the quadratic eigenvalue problem 
with the new eigenvalue parameter c in ([2]), 

9,9 x 40 2 4/1 1 \ , 

&o(fc 2 + L)Ur = ~^cU r + -{^---2 ;J cPU r . (14) 

Figure [5] shows numerical approximations of the first 
five curves of zero eigenvalues in the upper half of the 
(fii, ^2)-plane for fixed values of R\ = 1, R2 = 2, b a = 
0.4, and k = 1 and their mirror reflections in the lower 
half plane. The dotted curves show the diagonal line 
ill = ^2, the Rayleigh line Q\R\ — Q2R2, as wen as 
the axes Oi = and f^2 = 0. It is clear that each curve 
approaches the diagonal line = SI 2 for large values of 
Sli. 2. When bo becomes small, they approach closely to 
the line Sli = f^- 



The above conclusions also follow from rigorous analy- 
sis of the quadratic eigenvalue problem (|14[) . In the limit 
SI2 — > 00, we can set A = f^c as a new eigenvalue and 
treat the last term in (fT4"]) as a small bounded perturba- 
tion to the unbounded operator. In the limit 60 ~> 0, we 
set c = 6gA and again treat the last term in (|14p as a small 
perturbation. In both cases, eigenvalues A approach to 
the first eigenvalues of the positive unbounded operator 
r 2 (k 2 + L). We note, however, that this approximation 
is not uniform for all bifurcation curves and only apply 
to the finitely many bifurcation curves. 

To summarize, we gave mathematically rigorous proofs 
about distributions and bifurcations of eigenvalues of lin- 
earized operators associated with an ideal hydromagnetic 
CT-flow that lay a firm basis for identification of unstable 
modes in MRI experiments with real dissipative liquids. 
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